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Abstract 

Crackling noise is observed in many disordered non-equilibrium 
systems in response to slowly changing external conditions. Exam- 
ples range from Barkhausen noise in magnets to acoustic emission in 
martensites to earthquakes. Using the non-equilibrium random field 
Ising model, we derive universal scaling predictions for the dependence 
of the associated power spectra on the disorder and field sweep rate, 
near an underlying disorder-induced non-equilibrium critical point. 
Our theory applies to certain systems in which the crackling noise 
results from avalanche-like response to a (slowly) increasing external 
driving force, and is characterized by a broad power law scaling regime 
of the power spectra. We compute the critical exponents and discuss 
the relevance of the results to experiments. 



1 



1 Introduction 



Many non-equilibrium physical systems ranging from disordered ferromag- 
nets to superconducting vortices Q, to martensitic shape-memory alloys 
IfJ. [| [|, |], to earthquakes, exhibit crackling noise in response to smoothly 
varying external conditions (driving force). While it is not at all obvious 
that these systems should show similar properties, it is observed that these 
and other examples in nature exhibit power law scaling of the noise statistics 
over many decades. The exponents observed in these systems fall into distinct 
universality classes: Barkhausen noise observed in disordered ferromagnets 
exhibits power law scaling in the power spectrum (PS) 

P(u) ~ uj- 1/auz , (1) 

for large frequencies lu, with a universal exponent 1/avz that does not de- 
pend on the microscopic details, but only on a few basic properties such as 
symmetries, dimensions, range of interaction, etc. Since universal power law 
characteristics are often associated with systems at or near an underlying 
critical point, models with critical points have been suggested to understand 
the origin of the observed universality. 

Barkhausen noise, the characteristic crackling noise associated with the 
motion of magnetic domain walls in a ferromagnet as the external magnetic 
field is slowly varied, has enjoyed a significant amount of experimental and 
theoretical attention in the past decade because it is a particularly sim- 
ple and experimentally readily accessible example of systems with crackling 
noise. Models for Barkhausen noise in soft ferromagnets usually model the 
motion of one (or few) domain walls, while other models, presumably for 
Barkhausen noise in hard (strongly disordered) magnets model the collective 
behavior of many interacting domain walls. While the former is marked by 
long range ferromagnetic interactions which cause the system to "self tune" 
to the domain wall depinning transition |?[] the latter is not self tuned and 
requires a different mechanism to explain scaling. This mechanism comes in 
the form of a disorder induced critical point found in the zero temperature, 
random field Ising model (RFIM) || || studied far from thermal equilibrium 
[ TO | . In this paper we investigate the power spectra (PS) of "Barkhausen 



noise" in the zero temperature non-equilibrium RFIM and subsequent scal- 
ing relations near the disorder induced critical point. Many of the methods 
will be applicable to other systems with crackling noise 0, arising from the 
system's proximity to an underlying non-equilibrium critical point. Within 
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the RFIM the pulses of the Barkhausen signal are understood as collective 
events (dubbed avalanches) in which many magnetic domains (spins) flip in 
time as the external field is slowly ramped up or down. The properties of the 
Barkhausen signal and corresponding scaling relations are then equivalent to 
the statistical properties of the avalanches which are obtained by virtue of 
the system's proximity to a critical point. 

The association of an avalanche in the RFIM with a pulse in the Bark- 
hausen signal presents some practical complications in many experimen- 
tal systems: background noise and finite field sweep rates lead to merging 
avalanches and ambiguous definitions of pulses. It therefore seems best to use 
spectral tools [11, 12] to analyze the pulse train as a whole, rather than, say, 
trying to extract pulse time, or pulse size distributions. The Power Spectra 
(PS) that will be studied in this paper is defined as the absolute amplitude 
square of the Fourier Transform of the voltage time series V(t). The voltage 
V(t) obtained in Barkhausen noise measurements is the voltage induced in 
a pickup coil as a function of time t. The signal V(t) is proportional to 
the change in magnetization during the microscopic time interval from t to 
t + At, where At is the time it takes a single spin to flip. An alternative 
and very useful approach to the study of Barkhausen noise at finite sweep 
rate is provided by the ABBM model |13| , which is similar to the mean field 
approximation of the RFIM. As we will see, for both the ABBM model and 
the RFIM the high frequency scaling behavior of the PS is independent of 
the sweep rate and scales with a universal exponent 1 / {avz) having a value 
of 2 within the ABBM model and becomes around 1.77 within the RFIM. 

The organization of the paper is as follows. In section |] we define the 
zero temperature, non-equilibrium RFIM and review previous results on the 
scaling behavior of the avalanche size distribution and magnetic hysteresis 
curves near the underlying disorder induced critical point for an adiabatically 
slowly increased external magnetic field. In section || we introduce the PS, 
and give scaling predictions for the PS as a function of the amount of disorder 
near the critical point, again for the adiabatic case. In section ^ we generalize 
the results on the PS to the case where the external field is swept at a 
finite rate and give results for the expected critical exponents obtained from 
Widom scaling collapses of the numerical simulation results. In section |5| an 
application of the PS to detect lack of causality in the order of spin flips 
during an avalanche is presented. Finally, in section |] we summarize the 
results obtained and discuss related questions to be addressed in the future. 



3 



2 Hysteresis Modeled with Ising Spins 



In the zero-temperature RFIM magnetic domains are represented by Ising 
spins (sj = ±1) on a hyper-cubic d- dimensional lattice. Spins interact fer- 
romagnetically with their nearest neighbors with a strength J, and are sub- 
jected to a homogeneous external field H(t). Structural disorder is included 
by a quenched random field hi at each lattice site, with a Gaussian distribu- 
tion 

where R parametrizes the strength of the disorder. Spins interact with an 
Ising Hamiltonian 

n = -j^2s iSj -^2(H + h i )s i (3) 

The notation () implies summation over nearest neighbor pairs only. The 
temperature is set to zero and each spin is aligned with its local effective 
field /if , defined by 

hf = J E Sj + H + hi , (4) 

je<i> 

where < i > denotes the nearest neighbors of site %. The system is studied 
far from thermal equilibrium, i.e. it is typically not in the ground state. 

A dynamics is introduced into the system as follows: Initially the external 
magnetic field is at — oo, so that all spins are pointing down (sj = —1 for all 
% = 1, N). The magnetic field is then slowly increased. The effective field 

rr rr 

hf at each site is computed. If hf at site % changes sign, the corresponding 
spin Si at that site is flipped. A spin flip may trigger neighboring spins to 
flip as well, thus leading to an avalanche of spin flips, which is the analogue 
of a Barkhausen pulse in real magnets. 

In this paper, the external field will be ramped up and down at a constant 
rate 

°-£- 

3 The Adiabatic Limit 

The zero sweep rate limit Q — > is the adiabatic limit. In this limit the 
external magnetic field is kept fixed during an avalanche. Only after the 
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avalanche has come to a halt, the external field is increased until it triggers the 
next avalanche. Combined analytical [|10j and numerical approaches [14] have 
shown that in this case there is a disorder driven dynamical phase transition 
at a critical disorder R c (R c = 2.16 in three dimensions, for the RFIM with 
nearest neighbor interactions and Gaussian disorder, in units of the exchange 
coupling J), separating a low disorder regime (R < R c ) characterized by 
hysteresis loops with a macroscopic jump (AM) in the magnetization M, 
from a high disorder regime (R > R c ), in which the hysteresis loops look 
smooth, see Fig. [3]. Here the magnetization is defined as M = (X)i=i S i)/N, 
where N = L d is the total number of spins in the system, L its linear size and 
d the dimension. The simulation results reported in this paper are exclusively 
for d = 3. 





l/o 


T 


(35 


V 


H c M c 


0.035(30) 


4.20(30) 


1.60(6) 


1.8(2) 


1.4(2) 


1.435(4) 


0.9(1) 



Table 1: Critical exponents and critical fields defined in the text, obtained 
from simulations of the adiabatic case in d = 3. The values of H c and M c 
are not universal. Values quoted correspond to a RFIM with n.n interactions 



and with a Gaussian distribution of random fields. From 15 . 



The jump AM in the magnetization for R < R c scales to zero as 

AM ~ (R c - Rf , (6) 

where /3 is a universal prediction of the model, table |I[ At the critical dis- 
order (R = R c ) each branch of the saturation hysteresis loop has one point 
(M C (H C )) where the slope diverges, dM/dH(H c ) — > oo: near that point the 
magnetization is described by a power law of the form 

\H-HA 1 ' 8 



\M(H)-M C 



(7) 



where 5 is another universal prediction for experiments (table [l]). 

The apparently smooth parts of the hysteresis loops really consist of many 
microscopic steps not resolved in Fig. [I]. These steps are the avalanches of 
spin flips, in analogy to Barkhausen pulses in real materials. For R > R c 
the distribution of avalanche sizes D(S, R, H) (which is proportional to the 
probability to observe an avalanche of S spin flips at disorder R and external 
field H in response to a small increase in H) scales as 

D mt (S, R, H) ~ S- T V + (S°r, h/r? 5 ) , (8) 
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Figure 1: Right branches of hysteresis loops (for increasing H) for high dis- 
order R = 3, close to critical disorder R = 2.25 and low disorder R = 1.6. 
Results are for L = 200 at = 10" 6 . 

with r = (R — R c )/R c , h = (H — H c )/H c , and T> + is a universal scaling 
function [fTEl]. There are similar scaling forms for the avalanche correlation 



function, the cluster correlation function, and other quantities [|Iq| . The 
correlation length £ scales as the diameter of the largest avalanche of the 
power law scaling regime of D(S, R, H): 

$,{r,h)~ r - v E{h/rP 5 ) , (9) 

where S is a universal scaling function ||16||. Values for the exponents are given 



in table [l|. Interestingly, numerical simulations indicate that the "critical 
region" is remarkably large: almost 3 decades of power law scaling in the 
avalanche size distribution remain when measured at a disorder R that is 40% 
away from the critical point. At 2% away one extrapolates seven decades of 
scaling [TfJ. This may explain why in many experiments it does not even 
seem to be necessary to tune the disorder to see the critical power law scaling 
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over several decades: the samples used may just fall into this large critical 
region. 



3.1 Power Spectra 

Each avalanche, or Barkhausen pulse, has an internal structure reflecting 
its evolution with time, see Fig. |2|. We denote with V(t) either the voltage 
measured in a Barkhausen noise experiment, or the number of spins triggered 
in our model during the same microscopic time interval of length At, which 
is the time it takes a single spin to flip. In our model we set At = 1. The PS 
is defined as the amplitude square of the Fourier Transform of the voltage 

to+T 

P(u) = | e iut V{t)\ 2 . (10) 

t=t 

There is some freedom in choosing the initial time to an d the total duration 
of the transform T. We will provide concrete criteria to fix these parameters. 

In the strict adiabatic limit (Q = 0), avalanches are separated in time 
by an infinite time, and therefore, only the PS of an isolated avalanche may 
be rigorously defined. The PS spectra may be defined from the following 
limiting process: The frequency Qa is defined as the largest sweep-rate such 
that no two avalanches overlap in time. The frequency Qa will be finite in a 
finite system and will approach zero with increasing volume. The adiabatic 
limit can then be defined as the infinite volume result of the PS as computed 
with the sweep rate Qa- A general scaling form for the PS in the adiabatic 



case may be derived from generalizing the result of Kuntz and Sethna (T? 
to the present situation 

P(u, h, r, L) ~ uo- llrjvz T{uo- l,{vz) r, h/r? 6 , Lr v )L d \AM\ , (11) 



where AM = M(t + T) - M(t ). It was found in Ij, that at criticality, 
for r < 2 (as in the RFIM at critical field and disorder), P(tu) ~ i J j^ 1 l avz for 
large u. On the other hand, for r > 2, at the critical point, for large u, one 
finds P(uj) ~ uj( T - 3 )/ auz , The scaling form T is therefore more conveniently 



split into two contributions |17j so that the general form for the PS becomes 
P(u,h,r,L) ~ uo- llavz Ti{uo- ll[vz) r,h/r ps ,Lr v )L d \AM\ (12) 

where T\ and JF 2 are universal scaling functions, which become finite for large 
frequencies. 
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4 Finite Field Sweep rate 



The above scaling forms were all obtained for zero field sweep rate. In the 
following section we take into account that in experiments typically the ex- 
ternal magnetic field is ramped up and down with a constant finite sweep 
rate Q = dH/dt > 0. In order to compare sweep rates at different volume 
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Figure 2: Time series for the noise signal in a ferromagnetic system. This 
corresponds to a time series for sweep rates of increasing magnitude from A 
to D. 

sizes, it is convenient to define 



n = 



(13) 



where N is the number of spins in the actual system and v defines a size 
independent sweep rate. In contrast to the adiabatic case ( Q — > ) where 
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the external magnetic field is kept fixed during an avalanche, at finite field 
sweep rate it can trigger new avalanches before currently running ones have 
petered out, see Fig. 0. 

Interestingly, the Barkhausen noise time series at relatively fast sweep 
rates Q (as in case D in Fig. |^, where basically no separate avalanches can be 
identified since most individual avalanches overlap in time) has the same PS 
as the low driving sweep rate cases A, B or C, at least for high frequencies 
u, and therefore results in the same high frequency noise scaling exponents. 
This highlights the usefulness of the power spectra: it allows to character- 
ize Barkhausen noise in situations where individual avalanches cannot be 
observed separately. 

To anticipate the effect of finite sweep rate we note the general form of 



the power spectra of an ensemble of independent pulses given in [18[ where 
it is assumed that there are no size correlations between pulses: 



PtotaM = a(H)Q 



' , t „ , u,2 ( In° D(AT)e iuAT dAT V 

(14) 

where () is the average over all the individual pulses, $(u;) is the Fourier 
transform of an individual pulse and D(AT) is the distribution of time inter- 
vals between independent, successive nucleation events (i.e. spin flips that 
in the adiabatic limit would be triggered only by an increase in the external 
magnetic field rather than some neighboring spin flips). a(H) is the number 
of nucleation events per unit field increase, yielding a(H)QT nucleations in 
a time span T. 

As shown in Fig. |3], for the RFIM with only n.n. interactions D(AT) is 
Poissonian, leaving the second term in brackets zero. Hence the resulting 
PS is simply proportional to the sum of the individual PS. This holds as 
long as any simultaneously propagating collective events do not overlap in 
space (otherwise the pulses can no longer be considered independent |p9f| ). 
Conveniently, this implies that the PS can be used to determine many scaling 
properties of the simultaneously propagating adiabatic avalanches at high 
sweep rates. 

It is important to note here that the introduction of long range (anti- 
ferromagnetic) dipole-dipole interactions, as present in soft magnetic ma- 
terials, would lead to distinctly non-Poissonian D(AT) which results in a 
nontrivial sweep rate dependence in P tot {uj). We leave a detailed study of 
this dependence for the future. The results of this paper are expected to 
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Evidence of statistical independance of nucleations 



Without Anti-Ferromagnetic Interaction 

fit — 




50 100 

AT separating nucleations (units 1/L ) 

Figure 3: Plot of the distribution of time intervals AT separating successive 
nucleation events. The continuous line is a fit to a Poissonian distribution. 
For comparison, a schematic figure for an anti-ferromagnetic interaction is 
shown. 

apply to hard magnets and other systems, in which long range dipolar fields 
are negligible due to sample geometry and other factors. 

4.0.1 The determination of H c 

Before testing the above scaling relations Eq. [12], and Eq. [14], we first need a 
criteria to pick up a sensible initial time t and total length T of the Fourier 
transformed Barkhausen train to compute the PS (see Eq. [H|). In order to 
apply Eq. ^ at the critical field (h = 0) we choose the Fourier transformed 
interval such that 

H c e[H(t ),H(t + T)] . (15) 

Furthermore, in order to have sensible statistics, one would like to include 
as many spins flips as possible in the time series, while keeping the deviation 
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from H c small enough so that Eq. |T2| with h = is still applicable. We 
typically choose the intervals such that h = (H — H c )/H c < 0.05. This way, 
it is found that a reasonable fraction of the spins are included in the time 
series and the range of recorded events spans all sizes. 



4.0.2 How fast is fast? 

We now compute the critical sweep-rate Q c , above which one observes changes 
of the PS due to spatially overlapping simultaneous events. Since the nucle- 
ations of avalanches are random in time (i.e. Poisson distributed) Eq. [14] 
implies that the adiabatic scaling result for P(uj) will hold as long as the 
spin flip order is maintained for the adiabatic avalanches. That is, as long 
as we are superimposing the adiabatic avalanches to construct the time se- 
ries we are assured of obtaining a power spectra that is the superposition 
of the PS of the individual pulses. The flip order begins to change when, 
during the propagation of an avalanche, new avalanches are triggered within 
the same volume that the initial avalanche would have eventually covered. 
The increased sweep-rate leads to "parallel processing" and the avalanche 
that would have taken T time steps to complete now takes T' such that 
T' < T for all Q > O c (T). The field sweep rate at which the expected num- 
ber of spins flipped in avalanches triggered within the volume of the primary 
avalanche is some finite fraction of the primary avalanche size is given for an 



avalanche of duration T as [IS 



n c (T) = S-^— (16) 

Since the power spectra is a superposition of all the individual avalanches 
there will not be any change in the PS as long as 

n<n c = (17) 

where T* is the largest duration observed which scales as T* ~ £ z . Since 
£ ~ L if r~ u > L or £ ~ r~ u if r~ u < L the critical field sweep rate has two 
scaling forms 

Vl c ~ L d f {T ~ 2 ~ auz) for £ >> L (18) 
Q c ^ £d f (T-2-<rvz) for ^ << L ( 1Q ) 
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We note here that if r — 2 — ovz > — 1, as is the case in the RFIM, then 
even if the avalanches are compact, i.e. df = d, the onset of sweep rate 
effects happens at a field sweep rate much faster than the largest sweep rate 
at which one can still measure separate pulses (since Q is measured in units 
of 1/L d ). Fig. § confirms the applicability of the adiabatic result to a large 
range of external sweep rates Q < Q c . 

Furthermore even for f2 > f2 c we can derive a frequency tu a (Q) = fi£ a! ^ 2 ~' r ) 
so that for u > u> a (Q) the power spectra are again described by the adiabatic 



results [pL9|] . Below this frequency the PS follows a cutoff function as a result 
of the sweep rate decreasing the duration of the larger avalanches (see Fig. |]). 
For uj « Q one finds the asymptotic u — > result 

lim P(uj) = L 2d \M{t + T) - M{T)\ 2 = L 2d \AM\ 2 . (20) 

Recalling the scaling relation JTTJ 

AM ~ r^Mihr' 135 , Lr v ) , (21) 
near the critical point, in the adiabatic case (Q — > first) 

lim P(u) = L 2{d - p/u) for r -> , h -> 0. (22) 



Let us recall that while in fCT] , the limit of small a; is taken after the adiabatic 



Q — > limit has been performed, in our simulations we will be working at 
finite sweep rates and therefore, our PS will always satisfy Eq. ^0] correctly 
(with non-zero AM). It will be shown, however, that in d = 3 dimensions, 
(3/v will be very close to zero and therefore, even the adiabatic limit scales 
very approximately as L 6 . 

4.0.3 Critical exponents 

The PS is now obtained from numerical simulations and the scaling predic- 
tion Eq. IT^I in section RTTJ is used to extract critical exponents. The results 



presented involve averages over several, usually about 16, disorder realiza- 
tions. 

Besides the external driving frequencies, there are two additional natu- 
ral frequencies in the system defining the characteristic times of the largest 
events in the system: 

■* s w- ' (23) 
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Figure 4: Plot of the PS (for / = j-) as a function of Q at R c for L = 
100. Result includes Q = 10" 9 and = 10~ 6 (thin solid lines), £1 = 10" 4 
(dot-dashed line) and Q = 10~ 3 (wide solid line). The fit P(f) ~ f-Vi™*) 
corresponds to VI — 10 -9 and has been separated for proper visualization. 

where £(r) is the correlation length, and 

= ^ , (24) 

where L is the finite linear size of the system. 

For future convenience, we introduce a simple continuous function uj c 
that asymptotically approaches the values of u T and Ul, as £/L — > and 
£/L — > oo, respectively: 

_ l+g(r)«/L* 

= 77 r; = u; r + u;L. (25) 

Other definitions with the same assymptotics can also be used. 

In Fig. the PS for H from a small interval around H c is shown as 
a function of the disorder strength for fixed system size. As the disorder 
approaches R c , the system size L becomes the characteristic length scale 
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Figure 5: Power Spectra for (/ = at different disorder values R = 2.30, 
R = 2.50 R = 2.60, R = 2.70. The thin line corresponds to R = 2.20 < R c , 
and the dashed line is the fit to R = 2.30 (it has been shifted for proper 
visualization). Results correspond to L — 200, Q = 10~ 6 , and if in a small 
interval around H c with (H — H c )/H c < 0.05. 




Figure 6: Power Spectra (for / = 
to R = 2.7 at L = 200, SI = 10" 6 



2n' 



collapse for disorder values R = 2.3 



according to Eq. |27|. The value of A is a 
non-universal parameter dependent on the size of the system L. (A ~ L z ). 
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of the system implying that uo^ takes over u r as the characteristic cut-off 
frequency, and one observes a very large region of power law scaling (about 
four decades for L = 200). For larger disorder R > R c , it becomes apparent 
that oj c ~ oj r ~ r uz in agreement with the behavior suggested in Eq. 0. 
At large frequencies expected, with 

\/{<rvz) = 1.77(4) , (26) 

where the error bar merely indicates the dispersion around the mean value 
from the different fits to Eq. [12| corresponding to the different disorder pa- 
rameters and system sizes considered. 

More information can be obtained from a scaling collapse of the PS for 
different disorders. Setting h = in the adiabatic scaling form Eq. |T2| we can 
make an Ansatz for the specific form of the scaling function JF to arrive at 
the scaling form (assuming r < 2, as discussed) for j < 1: 

P(u>) = (— )-sQ(- )\AM\L d+l ' av . (27) 

The factor of | AM\L d is a normalization factor resulting from our specific 
definition of P{uj). Its particular form in the adiabatic case has been given in 
Eq. ^2[ The asymptotic behavior of the function Q is given by simple power 
laws (with x = ^-): 

«*) = {*£ • x> >] , (28) 

where the scaling for x » 1 implies the expected behavior ijj~ x l^ avz ) at large 
frequencies. Plugging Eq. EB! and Eq. ESI into Eq. H7\, for a; << 1, one obtains: 




P( W )„ + ,co«u c , (29) 



i.e. the PS becomes independent of the frequency. This is a consequence 
of the fact that the system is above criticality, there is a minimum cut-off 
frequency u r , which corresponds to the finite duration of the largest avalanche 
in the system, with no characteristic events having larger durations. This 
result is reflected by the low frequency plateau at large disorder in Fig. |[ 

The collapse obtained using the scaling form Eq. |2?| is plotted in Fig. ||. 
For the particular system size shown (L = 200), the collapse extends up to 



15 



2.5e+05 



2e+05 



^ 1.5e+05 



1e+05 



50000 



2.2 



2.3 



2.4 



vz=2.44 

(2-x)l(avz)=0.68 
R=2.25 




2.5 

R 



2.6 



2.7 



2.8 



Figure 7: The coefficient C(f r , fi) = 7^17^ , see Eq. as a function of dis- 
order i? for L = 200 and f2 = 10~ 6 (grey dots). The dashed line corresponds 
to a fit according to Eq. RTL 
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Figure 8: Finite size scaling collapse of Power Spectra (for / = at R — 
R c for system sizes L = 50 to L = 200, (L=50,75,100,150,200), 7r collapsed 
according to Eq. |34|. 
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about five decades and it only fails for very low frequencies that are of order 
or smaller than the (small) nonzero external sweep rate Q in the simulations. 
From the collapse one obtains the critical exponents 

— = 2.38(10) , = 2.41(10) , = 0.61(9) , (30) 

vz avz avz 

with error bars estimated as previously. 

It is found that L d AM is numerically very close to L d , and from this and 
the scaling form Eq. |2T] one derives that (3 is close to zero within the accuracy 
of our calculation. 

Additional information may be obtained by further exploiting the large 
frequency behavior in Eq. ^7], to » uj c . We can write 



2-T 

avz 



C(u r ,u L )=A(L) I— -j- , (31) 



1 + 



One can fit the coefficient C for different disorder values (see Fig. [7]), where 
L and therefore A(L) is a constant). The fit is excellent at the two ends 
of the curve and becomes slightly inaccurate on the region £/L ~ 1, which 
is a consequence of the approximation made in selecting the interpolation 
function for u c given in Eq. Since the exponents are computed from the 
two ends, one obtains 

— = 2.44(10) , = 0.68(8) , (32) 

vz avz 

which provides a cross check of the validity of the results obtained in the 
collapse. It is also interesting to compute the correlation length 

m = 9 - 68 v^r) ■ (33) 

The large exponent v ~ 1-4(1) is partly responsible for the large scaling 
region of the non-equilibrium RFIM [110 . 



4.0.4 Finite Size Scaling 

The analysis of the PS as a function of the linear system size L at criticality 
r = h = is also very interesting. From Eq. O one derives the adiabatic 
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scaling form 



P(uj)L 1/{au)+d \AM\ = P(uj)L d - p/v+l/{(7y) = g(ioL z ) , (34) 

with the universal scaling function G(x) ~ 1 for x — > and Q(x) ~ x~ x '^ vz ^ 
for x — > oo. Fig. |8] shows that the data collapse very well to the assumed 
relation for all the included frequencies (about five decades of scaling for 
L = 200). The resulting exponents are 

^ = 1.63(10) , l/(avz) = 1.75(6) , 9/2 + l/av = 3.0(5) (35) 



where the exponent relation 9 = d — (3/v — 1/ (erz/) [10 has been used and the 
error bars indicate the dispersion from the values obtained as a consequence 
of disorder fluctuations. 



4.1 Below R c 
4.1.1 General Forms 

The hysteresis loops of the RFIM below the critical disorder R c have a macro- 
scopic jump in the magnetization at a critical field H C (R) where many spins 
flip almost simultaneously in an event that extends over a finite fraction of 
the system. The jump grows with decreasing disorder. Studying the critical 
behavior of the RFIM below R c gets complicated by the fact that even if 
the correlation length becomes smaller than the system size, there is always 
an event which is sensitive to the system size, and in previous studies, large 
finite size effects have made the analysis at R < R c difficult. The PS below 
but near R c for uo > Q can be collapsed using the scaling form 

p( u ) = 0( UJ/UJL I )\AM\L 1 K av ') +d . (36) 

1 + UJ r /UJ 

with 4>(x) ~ x~ x l{ GVZ } for x » 1. The latter result states that short-time 
propagation of the avalanches proceeds in the same way, whether above or 
below (but close to) the critical disorder. At zero disorder, the PS is a 
constant: 

hm P(u) = L 2d . (37) 

According to Eq. |20], this also scales as the limit of the PS for R > for 
uj « Q. The total duration T&m(L, R) of the jump in the magnetization 
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for R < R c decreases for decreasing R. This leads to an increasing cut-off 
frequency 

^-Tzkwr (38) 

with P(u) ~ for to > u(r, L) and P{uo) -> CL^ 1 ^ |AM| with C 

a constant, for << uo « u{r, L) and P{oS) — > L 2d for c<j << f2. 

Note that the scaling form Eq. |36] is really a special case of the general 
form Eq. [12] for r < 2 derived above R c with different scaling functions T\ 
and JF 2 . 

Far below R c for very large systems we expect Eq. |36] to be replaced by a 
scaling form that involves the critical exponents of single interface depinning 
HI, |J, some of which are very close to the exponents of the critical point 
studied in this paper (for example, — is the same within error bars) [|H]]. 

For the numerical analysis we note that the general considerations dis- 
cussed in subsections [4.0. 1| and |4.0.2j concerning choosing the interval so that 
it includes H C (R), its width and the range of sweep rates, applies in the 
regime below R c as well. 



4.1.2 Critical Exponents 

The first issue we investigate is the large frequency limit, which is described 
by Eq. |36|. In Fig. | the results of the PS are plotted for different disorder 
parameters R < R c , for a low sweep rate Q = and in the adiabatic case 
(Q = 0) the fit to the large frequency behavior gives 

\j{pvz) = 1.79(5) , (39) 

with error bars defined as previously. 

Fig. |9| shows very clearly, particularly for the PS at R = 1 << R c , that for 
decreasing disorder P(uS) tends to a constant over larger and larger regimes 
uo < u(r, L), and consequently the large frequency region where the behavior 
Eq. |36] holds, shrinks. 

Fig. |TI] shows a collapse of the PS for different disorders according to 
Eq.(^) over about six decades of scaling for L = 200, with the critical 
exponents 

uz = 2.38(10) , l/{avz) = 1.77(5) . (40) 

In contrast to the analysis for R above R c , the analysis below R c does not 
involve the exponent r, related to the distribution of avalanche sizes. 
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Figure 9: The Power spectra (for / = for data below R c for L = 200 
and disorders between R = 2.2 and R = 1.6. R = 1.0 is plotted with a thin 
dotted line, for Q = 10~ 6 . All other curves are plotted for Q = 10~ 6 and 
Q = (adiabatic limit). The fit (straight dashed line) is slightly shifted for 
proper visualization. 
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Figure 10: Collapse of the Power spectra (for / = j-) for disorders R = 2.2 
to R = 1.6, for L = 200 and Q = 10 -6 . The data is collapsed according to 
Eq. |36], and f L = u L /27i and f r = u r /2ir. 
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More information below R c may be obtained from the region £(r) >> L 
where finite size effects dominate. Actually, the PS may be collapsed exactly 
as in Eq. |33|, and the exponents 

z = 1.63(10) , l/(avz) = 1.75(6) , (41) 

follow. Since the PS in this region becomes insensitive to the correlation 
length £, the results are identical to the ones above R c . 

5 The Interpretation of the Power Spectra 

The previous analysis has shown that the PS is sensitive to the causality 
relations among spins. This point can be made more explicit by considering 
a different dynamics than the causal dynamics implemented in this paper. 
Let us consider a standard metropolis dynamics, as if we were introducing 
temperature into the model. The algorithm is as follows: Spins are randomly 
selected^, and the standard Metropolis (or Glauber) acceptance/rejection 
check is used to determine whether the spin is to be flipped. Let us recall 
that since we are considering the temperature to be zero, the proposed move 
is accepted if the local energy = —hl"<Ji of the spin being tested is positive 
and rejected if it is negative. After N such attempts to flip randomly chosen 
spins are performed, the time is increased by one unit. 

A spin can only flip when it is selected, and therefore, if that particular 
spin is not selected by the random update, there will be a delay between the 
time when it is energetically favorable (which is the time when the spin would 
be flipped in the causal dynamics) and the time in which the spin actually 
flips. Spins flipping at earlier times as compared to the causal dynamics 
may also take place, since the random process may flip a spin seeding an 
avalanche, and within the same unit of time during the N attempts before 
updating the time by one unit, some of the neighbors to this seeding avalanche 
spin, which become now energetically unstable, may be flipped. Within the 
causal dynamics those spins would flip at later times. 

In Fig. |Tl| the time series generated from these two dynamics are shown 
and are extremely close. The PS, however, as shown in Fig. [L^, has a different 
high frequency behavior, and it is clearly sensitive to the different dynamics. 

1 since the disorder is itself random, a sequential update has very similar effect as a 
random one, with very similar PS 
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Figure 11: Voltage time series for causal and random updates for a sin£ 
disorder realization R = 2.3). 
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Figure 12: Power Spectra (for / = j-) of the previous time series. 
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Of course, time averages of the magnetization, etc.. will be insensitive to 



the different dynamics, as is clear from the time series Fig. |TT[ Observables 
that show some universality at large frequencies, as is the case for the PS, 
become sensitive to the different dynamics of spins, even if averaged over 
several disorder configurations. 

The previous remarks are of interest for extensions of the previous study 
to more general situations, in particular for considering temperature fluctua- 
tions when the system is out of equilibrium. Temperature would be naturally 
introduced into the system by considering Glauber dynamics with random 



updates, as done in other simulations with Ising spins f2T| . Of course any 
algorithm to study temperature effects should meet the necessary require- 
ments: 

• Converge to equilibrium for very slow sweep rates at finite temperature. 

• Reduce to the causal dynamics for small temperatures at finite sweep 
rate. 

It follows from our results that the second condition is clearly violated for 
random updates (or any other update uncorrelated with the dynamics of the 
system, like a sequential update). The update must be chosen correlated with 
the dynamics, thus reflecting the causal nature of the events (the avalanches). 



6 Conclusions 

6.1 Summary of results 

In this paper we have seen that the PS is an extremely valuable observable 
for the study of collective noise in physical situations where the avalanche 
picture would no longer be applicable. 

Besides its obvious theoretical interest, the PS becomes a powerful prac- 
tical tool to compute critical quantities. We have computed the following 
exponents summarized in table 0. 

The exponent predictions lie well within the error bars of most Bark- 
hausen noise exponents from experiments quoted in the literature, ||. Be- 
sides, we are able to obtain accurate exponents for the region below R c . 

The exponent quotes have been deliberately computed at relatively fast 
sweep rates, see section |] for a detailed discussion, thus providing conclusive 
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Above R r 



Below R r 



Quantity 


PS 


RFIM 


PS 


RFIM 


l/{avz) 


1.77(3) 


1.75(9) 


1.79(6) 




z 


1.63(10) 


1.7(2)* 


1.63(10) 




vz 


2.38(10) 


2.39(40)* 


2.38(10) 




3-t 
uvz 


2.41(10) 


2.45(8)* 






•2~t 

CTUZ 


0.61(9) 


0.70(6)* 







Table 2: Exponents computed and the ones known 

evidence of the irrelevance of low enough finite field sweep rate for the PS. 
Finally, it has been shown that the PS is sensitive to the particular dynamics 
of the spins. 



6.2 Experiments and Outlook 

In Fig. |i~3|our results are compared to the experiments from ||21|| , correspond- 
ing to a FeSi sample. Our result for the large frequency dependence of the 
PS looks more plausible than the exponent — = 2, suggested by the ABBM 
model ]n|. In [p^] , the exponent — = 2 is calculated reporting values in 
the range — = 1.5 — 2, in agreement with our results, but with a too large 
uncertainty for a precise comparison. Other experimental determinations in 
amorphous ribbons of Fe^Co2\Bvo also seem to favor a value close for — 
close to 1.78 [p3fl . For this particular data we expect that the frequency de- 
pendant bump in the PS is due to two things: (1) the presence of long range 
anti-ferromagnetic (LRAF) interactions which (change the nucleation dis- 
tribution from Poisson to something more complicated, see section |], and) 
affect the low frequency of the power spectra; (2) the fact that the larger 
avalanches take less time due to "parallel processing" discussed above also 
diminishes the adiabatic PS in the low frequency regime. While the full de- 
tails of these effects were not addressed in our current study we are presently 
investigating effects of LRAF interaction in the RFIM. 

In experiments on soft magnetic materials, Barkhausen noise seems to 
be due to the propagation of a single domain wall, and nucleation of new 
domains are rare events 0. Far below R c , our model also exhibits domain 
wall propagation during the macroscopic jump AM in the hysteresis loop. 
For small disorder R, the number of small nucleating events are expected to 
be negligible, and in this case the PS for the jump at high frequencies should 
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Figure 13: Experimental results from |2l). The dotted line is our prediction 
1.78. the dashed line is the result of 2 from the ABBM model []I3|. 



match those of single interface propagation in the absence of LRAF [|| [l7j 
(In some experimental systems the LRAF can be suppressed, for example by 
applying stress to the sample fl23"|). It is known that the critical exponents for 
single interface propagation belong to a different universality class, although 
the critical exponents characterizing the large frequency limit of the PS turn 
out to be numerically very close |], |7J. Therefore, a very detailed study 
is needed to discriminate between these two regimes (low and near critical 
disorder) in which the two models apply to the PS. 

The main problem for a detailed comparison of the predictions for the 
disorder induced critical point presented in this paper with the experiments, 
is the lack of control over the structural disorder in experimental systems. 
Recent experiments on Co/CoO |24j have shown that in these systems it is 
possible to control the amount of structural disorder experimentally. If Bark- 
hausen noise is measured in these disorder controlled samples, one could use 
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the scaling forms derived in this paper and get a very detailed understanding 
of the universal aspects of hysteresis, including detailed critical exponents 
near the critical disorder. 

Other experimental systems for which such analysis could be extended 
include acoustic emission in martensitic materials in response to stress or 
temperature ramping. In fact, the Barcelona group 0, ^|, [|, [5| has shown 
that crackling noise in martensites displays critical scaling in avalanche size 
distribution and the PS similar to the ones described in this paper for Bark- 
hausen noise, although with different exponents and therefore belonging to 
a different universality class. 

The analysis presented in this paper certainly allows several important 
applications which we have briefly discussed. First of all, the irrelevance of 
the sweep rate has very important practical consequences for future calcula- 
tions since by using a faster sweep rate one can compute the same quantities 
with less effort. This is an important benefit in situations, as for example the 
inclusion of temperature, where slow sweep rates are computationally very 
demanding. 

Another important aspect is how long range dipole-dipole interactions 
may affect the results, a subject briefly discussed in Fig. [3[ We hope to 
report more on that in the near future. 
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